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Abstract. Here the results of calculations of pulses 
in Cherenkov light detectors for the Yakutsk array 
are presented. As long as the Vavilov-Cherenkov light 
is used to calibrate signals in scintillation detectors at 
the Yakutsk array it is vital for these measurements 
to be precise. The vaUdity of measurements of the 
signal Q(400) used as the estimation parameter at the 
Yakutsk array has been confirmed. Our calculations 
show that the width of time pulses increases from 
nearly 50 ns at a distance of 100 m from the shower 
axis up to 700 ns at 1000 m. 
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I. Introduction 

It is a known problem of calibration of energy esti- 
mations of extensive air showers (EAS). At the Yakutsk 
array the Vavilov-Cherenkov light is used to calibrate 
signals in scintillation detectors Ul. Although our calcu- 
lations show that for the same flux of Vavilov-Cherenkov 
light there should be much higher signals in scintillation 
detectors than it is measured in the experiment ||2l- 

First, this contradiction may have different explana- 
tions. The hadron interaction model may not correctly 
describe the real processes at high energies. This ques- 
tion will be probed at LHC once it is operational. Alter- 
natively, we may have a different chemical composition 
of primary cosmic rays - all our calculations were under- 
taken in the assumption of light chemical composition 
at high energies. Although recent experiments claim that 
the primary particles at those energies are protons [31, 
m (with exclusion of |5|). 

Each detector array encounters the problem of full 
signal registration. Each detector after being triggered 
by some event (a particle hits the detector, the counting 
rate exceeds some level, a signal comes from another 
detector, etc.) collects the signal for some period of 
time, a so-called time gate. Time gates should be wide 
enough to collect possibly all particles coming to the 
detector On the other hand, since we always have a 
background from low energy cosmic rays or the local 
sources (radioactivity, light pollution, etc.), the time 
gates should not be too wide in order to keep the signal 
to noise ratio high enough. So the time gates should 
be of about the disk thickness. At most detector arrays 
such as Haverah Park, Volcano Ranch, Yakutsk, the time 
gates were set to about 2 /is. 

A. Watson in ||6l suggested that the time gates at 
the Yakutsk array are too narrow to detect all particles 
and thus the signals are underestimated, that leads to 
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Fig. 1: Pulses of Vavilov-Cherenkov light at 400 m dis- 
tance from the shower axis. Showers induced by primary 
gamma with energy 10 GeV starting from different depth 
(solid curve 1 - 350 g/crn^, dashed curve 2 - 550 g/cnP, 
dotted curve 3 - 750 g/cm^) 
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Fig. 2: Pulses of Vavilov-Cherenkov light at 400 m 
distance from the shower axis. Showers induced by 
primary electron with energy 316 MeV starting from 
different depth (solid curve 1 - 350 g/cm^, dashed curve 
2 - 550 g/cm^, dotted curve 3 - 750 g/cm^) 



an overestimation of energy of the primary particles. 
A rapid change in steepness of the measured lateral 
distribution functions of signals and the fact that some 
of the showers detected at Haverah Park had signal 
width more than 2.2 //s ||6l are claimed to support this 
suggestion. 

In our previous paper fT\ we have studied the problem 
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Fig. 3: I'ulses of Vavilov-Cherenkov light at 400 m dis- 
tance from the shower axis. Showers induced by primary 
muon with energy 10 GeV starting from different depth 
(solid curve 1 - 350 g/cm^, dashed curve 2 - 550 g/cm^, 
dotted curve 3 - 750 g/cm^) 
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Fig. 4: Fraction of Vavilov-Cherenkov light collected 
within time r at the different distances from shower axis 
(doted curve - 100 m, dash-dotted curve - 400 m, dashed 
curve - 600 m, solid curve - 1000 m) 
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of form and width of pulses in scintillation detector In a 
few words though the signals in scintillation detectors at 
large distances from the shower axis (more than 1 km) 
exceed 2 /is (at 600 m signal width is well within 2 /is) 
and only part of the signal is measured the difference 
between the theory and the experiment is still too high 
to be explained solely by partial signal measurements. 
Therefore, additional analysis was required. In this paper 
we investigated another question - the form of the 
Cherenkov light pulse at the detector. 

II. Method of simulations 

As well as in the previous work in this study we 
used CORSIKA 6.500 |8| with QGSJET-II |jl-|U3| for 
high energy and Gheisha-2002d lfT4l for low energy 
calculations. Parameters of the atmosphere and the mag- 
netic field were set to fit the conditions of the Yakutsk 
array. The time delay of arrival of the Vavilov-Cherenkov 
photons was calculated as the time between its arrival at 
the observation level and the time of arrival of the very 
first particle at the observation level. All calculations 
were carried out in the framework of a multilevel hybrid 
scheme |15 |, i.e. we calculated the database of Vavilov- 
Cherenkov light pulses at different distances from the 
shower axis from low energy showers from primary elec- 
trons, gamma and muons (some examples are shown in 
Fig. 1-3) and using source functions (also calculated in 
our previous simulations |7|) we calculated the impulses 
for high energy showers. 

Source function describes the number of particles with 
energy in interval \E\E + dE\ born in the atmosphere 
at depth \x\x \- dx\. The database of the Vavilov- 
Cherenkov light pulses Ce,'f,fj,{T,r, E,x) describes the 
Vavilov-Cherenkov photon flux from a shower induced 
at depth x by electron or gamma or muon with the 
energy E at the distance from the shower axis r at 
moment of time t. Time starts with the arrival of the 



very first particle at the observation level. Having the 
analytical source functions of low energy particles from 
ultra high energy shower and databases of the Vavilov- 
Cherenkov light pulses from low energy showers one 
can calculate pulses from ultra high energy shower: 

Fir,T) = Jdx r SeiE,x)Ce{T,r,E,x)+ 

XO Ethr 

+ J dx J S^{E,x)Cj{T,r, E, x) + 

+ Jdx Y S^{E,x)C^,{T,r,E,x). 

Here Ethr - some reasonable maximum energy thresh- 
old, maximum energy of particles in source functions; 
Emin - some low energy threshold (naturally in case of 
the Vavilov-Cherenkov light for electrons and gammas it 
is 21 MeV, for muons - 4 GeV); Se,-y,^{E, x)) - source 
functions of electrons, gammas and muons respectively. 
But if the source functions were obtained numerically 
as a table of parameters (energy Ei, depth of generation 
Xi and weight Wi if thinning procedure was applied) of 
low energy particles generated in the shower than this 
integral converts to a sum: 

Fi'r,T) = J2wiC^eiT,Ei,Xi,r)+ 

i 

+ J2wjCy{T,Ej,Xj,r)+ 
j 

+ J2wkCfj.{T,Ek,Xk,r). 

k 

After all theese calculations for vertical showers from 
primary protons we obtained the following results. 

III. Results 

In Fig. 4 the fraction of Vavilov-Cherenkov light 
collected within time r at different distances from the 
shower axis is shown. As one can see the width of 
the Vavilov-Cherenkov light impulse at the distance of 
100 m from the shower axis is less than 50 ns, at 400 m 
- 200 ns, at 600 m - 400 ns and at 1000 m it is 
less than 1 /is. Therefore, the time gates at the Yakutsk 
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Fig. 5: The Vavilov-Cherenkov light forefront (squares) 
and backfront (circles) and their approximations 

array detectors are set wide enough to collect all of the 
Cherenkov photons from EAS, thus in the experiment 
the Vavilov-Cherenkov light is measured correctly. 

In addition we have calculated the form of the 
Vavilov-Cherenkov light forefront and backfront (see 
Fig. 5). As the forefront arrival time we took the 
moment when 5% of total signal is collected and the 
backfront was considered as the moment when 95% of 
all Cherenkov photons were collected. We found that the 
best fit is achieved by the power function a — a ■ 
(with a/ = 1.63 • IQ-^, bf= 1.71, 3.95 • 10"^ ^nd 
bb= 1 .43 for forefront and backfront respectively) rather 
than with a spherical front. 

IV. Conclusion 

The width of the Vavilov-Cherenkov light pulse does 
not exceed 1 /is at the distance of 1000 m from the 
shower axis. Therefore, at the Yakutsk array all the 



Cherenkov photons from EAS are collected and the 
time gates of the Cherenkov detectors may be set even 
lower. 

Cherenkov light forefront and backfront may be ap- 
proximated by a power function a = a ■ (with 
a/ = 1.63-10-^ bf= 1.71, at = 3.95-10-2 and bb= 1.43 
for forefront and backfront respectively) with a better 
accuracy than with a spherical front. 
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